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Abstract 


We present the direct-imaging discovery of a giant planet orbiting the young star AF Lep, a 1.2 Mẹ member of the 
24 +3 Myr 8 Pic moving group. AF Lep was observed as part of our ongoing high-contrast imaging program 
targeting stars with astrometric accelerations between Hipparcos and Gaia that indicate the presence of substellar 
companions. Keck/NIRC2 observations in L' with the vector vortex coronagraph reveal a point source, AF Lep b, 
at +340 mas, which exhibits orbital motion at the 6c level over the course of 13 months. A joint orbit fit yields 
precise constraints on the planet's dynamical mass of 3.270 M. Jup, Semimajor axis of 8.4*11 au, and eccentricity of 
034 077. AF Lep hosts a debris disk located at ^50 au, but it is unlikely to be sculpted by AF Lep b, implying 
there may be additional planets in the system at wider separations. The stellar inclination (i4, = 54*3!^) and orbital 
inclination (i, = 5075 are in good agreement, which is consistent with the system having spin-orbit alignment. 
AF Lep b is the lowest-mass imaged planet with a dynamical mass measurement and highlights the promise of 
using astrometric accelerations as a tool to find and characterize long-period planets. 


Unified Astronomy Thesaurus concepts: Extrasolar gaseous giant planets (509); Astrometric exoplanet detection 
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(2130); Direct imaging (387); Orbit determination (1175); Debris disks (363) 


1. Introduction 


A growing number of planets with masses spanning 1— 
13 Myyp and separations from ~5 to 100 au have been imaged 
over the past 15 yr (e.g., Marois et al. 2008, 2010; Lagrange 
et al. 2010; Rameau et al. 2013; Macintosh et al. 2015; Chauvin 
et al. 2017; Keppler et al. 2018). It remains unclear whether 
these long-period “super Jupiters” predominantly form through 
a bottom-up core accretion scenario (Pollack et al. 1996), 
possibly assisted by pebble accretion (Johansen & 
Lambrechts 2017); through a top-down gravitational instability 
route (Boss 1997); or even via dynamical capture in a dense 
cluster at an early age (Perets & Kouwenhoven 2012). Further 
complicating the picture, the star formation process is expected 
to produce objects with masses as low as the opacity limit of 
~3 My through cloud fragmentation (Bate 2009). Demo- 
graphic studies are hinting that imaged wide-separation planets 
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have a bottom-heavy mass distribution (Wagner et al. 2019), a 
semimajor axis distribution weighted toward smaller separa- 
tions (Nielsen et al. 2019), an eccentricity distribution peaked 
at low values in contrast to higher-mass substellar companions 
(Bowler et al. 2020; Nagpal et al. 2023), and preferentially 
aligned stellar obliquities (Bowler et al. 2023). These results 
resemble the demographics of closer-in planets found with 
radial velocities (RVs; e.g., Fulton et al. 2021) although the 
significance of this similarity is limited by the small number of 
imaged planets. 

An efficient strategy to identify long-period companions is to 
use the gravitational reflex motion they induce on their host 
stars. Beginning with the discovery of the brown dwarf 
companion HR 7672 B (Liu et al. 2002), RV trends have been 
repeatedly used to identify promising targets for follow-up 
high-contrast imaging (e.g., Crepp et al. 2014; Cheetham et al. 
2018; Bowler et al. 2021). More recently, astrometric 
accelerations using proper motion differences between Hippar- 
cos and Gaia have proven to be similarly advantageous as a 
strategy to efficiently search for and characterize substellar 
companions (e.g., De Rosa et al. 2019; Brandt et al. 2019; 
Franson et al. 2022; Hinkley et al. 2023; Currie et al. 2023). 

The Astrometric Accelerations as Dynamical Beacons 
program is a high-contrast imaging survey targeting a 
dynamically informed sample of young accelerating stars 
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Figure 1. Left: predicted mass as a function of separation for AF Lep b from its HGCA astrometric acceleration. This plot is generated following Franson et al. (2023). 
The blue curve shows the masses and separations most consistent with the astrometric acceleration. Due to orbital phase, companions can occupy separations and 
masses above but not below the curve. The dynamical mass and semimajor axis of AF Lep b is denoted by the gold star. The planet is consistent with causing the 
acceleration. Right: Keck/NIRC2 Imaging L’ imaging of AF Lep b in 2021 December (left) and 2023 February (right). Here, S/N plots are displayed to enhance true 
sources over residual speckles. Note that the two images do not share the same color map. We apply a Gaussian filter with a standard deviation of 1.5 pixels to average 
over pixel-to-pixel noise. AF Lep b is detected with a significance of 7.50 in the 2021 December epoch and 13.0c in the 2023 February epoch. 


identified in the Hipparcos-Gaia Catalog of Accelerations 
(HGCA; Brandt 2021). Our targets are selected for youth, 
proximity, and from individual predictions of the companion 
mass and separation required to produce the observed proper 
motion difference between Hipparcos and Gaia EDR3. The first 
discovery from this program was the low-mass brown dwarf 
HIP 21152 B (Franson et al. 2023), a T-dwarf companion in the 
Hyades cluster. 

Here we present the discovery of a giant planet orbiting the 
young star AF Lep'* based on observations with the NIRC2 
vector vortex coronagraph (VVC) at the Keck II telescope. In 
Section 2, we summarize the physical and astrometric proper- 
ties of AF Lep. Our NIRC2 observations and analysis of the AF 
Lep light curve from the Transiting Exoplanet Survey Satellite 
(TESS) are described in Section 3. In Section 4, the relative 
astrometry of AF Lep b is used to confirm common proper 
motion with its host star and carry out an orbit fit. We also 
compare the orbital inclination with the stellar inclination to 
assess consistency with spin-orbit alignment. A summary of 
our results can be found in Section 5. 


2. AF Lep: A Young Sun-like Star 


AF Leporis (—AF Lep, HD 35850, HIP 25486, HR 1817) is 
a bright (V = 6.30 + 0.01 mag; Høg et al. 2000) F8V star (Gray 
et al. 2006) in the @ Pic moving group (Zuckerman & 
Song 2004). It has a distance of 26.84 pc (Gaia Collaboration 
et al. 2021), an age of 24 + 3 Myr (Bell et al. 2015), and a mass 
of 1.20 + 0.06 Mo (Kervella et al. 2022). The properties of AF 
Lep are summarized in Table 1. An infrared excess based on 
WISE, Spitzer/MIPS, and Herschel/PACS photometry indi- 
cates that AF Lep has a debris disk, and spectral energy 
distribution (SED) modeling implies the disk radius is located 
at 462-9 au although its unresolved nature means it could 
conceivably lie anywhere from 30 to 70au (Pawellek et al. 
2021; Pearce et al. 2022). Using dynamical arguments, Pearce 
et al. (2022) inferred that a single planet with a minimum mass 
of 1.1 X: 0.2 Mj, at a maximum semimajor axis of 35 + 6 au 


14 AF Lep b was also independently and contemporaneously discovered by De 
Rosa et al. (2023) and Mesa et al. (2023), both using Very Large Telescope/ 
SPHERE. 


would truncate the disk's inner edge. The minimum planetary 
mass increases steeply as its semimajor axis decreases, so a 
more massive planet closer in could also interact with the disk 
in a similar fashion. In a multi-planet scenario, the required 
planet mass is much lower (0.09 + 0.03 Mjy,). AF Lep has 
been frequently targeted in direct-imaging exoplanet searches 
(Biller et al. 2013; Stone et al. 2018; Nielsen et al. 2019; 
Launhardt et al. 2020). These searches did not find any close 
companions. 

AF Lep is listed as an SB2 in the chromospherically active 
binary catalog complied by Eker et al. (2008). The initial 
double-lined spectroscopic binary classification was reported in 
Nordstróm et al. (2004), who found broad-line profiles with 
structure and variable central positions (B. Nordstróm 2022, 
private communication). This variability may be caused by a 
close binary but could also originate from strong imprints from 
starspot-related activity. Indeed, Zúñiga-Fernández et al. 
(2021) did not identify AF Lep as a single- or double-lined 
spectroscopic binary based on their recent high-resolution 
optical spectroscopy of this system. Additional spectroscopic 
monitoring would help clarify the binary nature of the host star. 
Because the original classification of AF Lep as an SB2 has not 
been validated with follow-up observations, for this work we 
treat AF Lep as a single active young star. 

AF Lep exhibits a significant? astrometric acceleration 
between Hipparcos and Gaia DR3 in the HGCA. The HGCA 
provides three proper motion measurements: the proper 
motions at the Hipparcos and Gaia epochs and a joint average 
proper motion from the difference in sky position between the 
two missions. The average? tangential acceleration between 
the Gaia proper motion and the joint proper motion is 
0.0191 + 0.0024 mas yr ?, which corresponds to a physical 
acceleration of 2.43 + 0.30 ms! yr! at the distance of AF 
Lep. The predicted mass of a companion causing this 
acceleration would fall in the planetary regime if located 
within 2230 au (see Figure 1). 


m x? — T1, which corresponds to 8.5c for 2 degrees of freedom. 


16 The average acceleration is computed across a baseline of 12.6 yr. Note, 
though, that the joint proper motion reflects the motion of AF Lep over the 
entire 25 yr between the two missions. 
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3. Observations 
3.1. Keck/NIRC2 Adaptive Optics Imaging 


We obtained high-contrast imaging of AF Lep on UT 2021 
December 21 and UT 2023 February 3 with the NIRC2 camera 
at W.M. Keck Observatory. Our observations were carried out 
in L'-band (3.426—4.126 um) with the VVC (Serabyn et al. 
2017) and natural guide star adaptive optics (Wizinowich 2013) 
with the Shack-Hartmann wavefront sensor. Imaging was 
taken in groups of 25 frames using the quadrant analysis of 
coronagraphic images for tip-tilt sensing (QACITS; Huby 
et al. 2015, 2017) algorithm to center the host star behind the 
VVC. Each subsequence includes an off-axis unsaturated 
image of the host star for flux calibration and sky background 
frames for both the science images and the flux calibration 
frame. The science exposures have integration times of 0.18 s 
with 120 coadds. We used a subarray of 512 x 512 pixels for 
shorter readout times. For the 2021 December epoch, our total 
integration time was 39.6 minutes, amounting to 35°8 of frame 
rotation. For the 2023 February epoch, our total integration 
time was 131.4 minutes, yielding 85°5 of frame rotation. The 
Differential Image Motion Monitor (DIMM) seeing during our 
2021 December observations averaged 075; the average DIMM 
seeing during our 2023 February epoch was 0/6. 

Our data reduction for both data sets begins with subtracting 
darks and flat-fielding. The L.A.Cosmic algorithm (van 
Dokkum 2001) is then applied to identify and remove cosmic 
rays and hot pixels. We apply the distortion solution from 
Service et al. (2016) to correct for geometric distortions in the 
optics of the imaging system. The sky background is modeled 
and subtracted from both the science and off-axis flux 
calibration frames using principal component analysis (PCA) 
with the Vortex Image Processing (VIP; Gomez Gonzalez 
et al. 2017) package. Frames are co-registered to the sub-pixel 
level through a cross-correlation approach described in Guizar- 
Sicairos et al. (2008) and implemented in scikit-image 
(van der Walt et al. 2014) and VIP. Absolute centering is 
performed by fitting a negative Gaussian to the vortex core. 

For carrying out PSF subtraction, we experimented with 
reductions from both pyKLIP (Wang et al. 2015) and VIP. 
Both packages use PCA (Amara & Quanz 2012; Soummer 
et al. 2012) to model and subtract the host-star PSF. For the 
2021 December data set, the highest signal-to-noise ratio (S/ 
N) is produced by a pyKLIP reduction with 20 annuli, 20 KL 
modes, and a value of one for the movement parameter. For the 
2023 February imaging, the highest S/N is produced by a VIP 
reduction with 25 principal components. For the reduction of 
the 2023 February data set, the PSF is subtracted in annuli, with 
frames only included in the reference PSFs for a given science 
frame if they have more than 1 FWHM of parallactic angle 
rotation for the annulus. For each sequence, we then compute 
S/N maps by measuring the flux in 0.5-FWHM-radius circular 
apertures. The noise level is estimated through the flux 
measured in nonoverlapping circular apertures at the same 
separation. Figure 1 shows the S/N maps for our two epochs. 
We detect AF Lep b at a S/N of 7.50 in the 2021 December 
imaging and 13.0c in the 2023 February imaging. 

To mitigate the introduction of systematics from the PSF- 
subtraction algorithm, we use the negative companion injection 
approach (e.g., Marois et al. 2008) to measure astrometry. A 
PSF template is generated by median-combining the off-axis 
flux calibration frames taken over each sequence. The inverse 
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of this is then injected in the preprocessed frames at the 
approximate position, separation, and contrast of the compa- 
nion, and the PSF-subtraction routine is performed. If the 
astrometry and photometry is well matched to the companion, 
its signal will be removed in the post-processed images. 
Otherwise, the companion's parameters are adjusted, the 
negative template is re-injected, and the PSF-subtraction 
routine is run again. 

The injected astrometry and photometry are first coarsely 
optimized using the AMOEBA downhill simplex algorithm 
(Nelder & Mead 1965). We then use the emcee affine- 
invariant Markov Chain Monte Carlo (MCMC) ensemble 
sampler (Foreman-Mackey et al. 2013) with 100 walkers and 
2 x 10* total steps to explore the astrometry and photometry 
parameter space. The sum of the residuals within a 4 FWHM- 
radius circular aperture of the source is used within the MCMC 
sampler to assess the match of the injected PSF to the 
companion. The first 30% of each chain is discarded as burn-in. 
We evaluate convergence through both visual inspection of the 
chains and by verifying that subsections of the chains yield 
consistent astrometry. For the 2021 December data set, the 
companion only appears at a high S/N in the PCA reductions 
when a parallactic angle cut is applied in constructing the 
reference PSFs for a given science image. We use 85 PCs for 
measuring the astrometry of the 2021 December epoch, which 
produces the highest S/N among the PCA reductions. For the 
2023 February imaging, the companion appears at a high S/N 
in reductions without a parallactic angle cut, so we do not 
include this cut in the astrometry measurement to increase 
computational efficiency. When a parallactic angle cut is not 
performed, the maximum S/N is produced by 16 PCs, so we 
adopt this for the measurement. Ultimately, the number of PCs 
and PSF-subtraction approach for the astrometry measurement 
and reduced images are each a balance between S/N, 
computational efficiency, and consistency in the approach for 
the two data sets. 

The MCMC chains consist of separation p, position angle 0, 
and a flux scaling factor for the PSF template. Following 
Franson et al. (2022), we combine the uncertainties on the 
astrometry from the MCMC run with the uncertainties in the 
distortion solution, north alignment, and plate scale. We also 
incorporate a 4.5 mas centering uncertainty to account for the 
average performance of the OACITS algorithm (Huby et al. 
2017). Our final astrometry for the 2021 December epoch is 
p = 338 + 11 mas and 0 = 62°8 + 1°3. Our astrometry for the 
2023 February epoch is p = 342 + 8 mas and 0 = 72°0 + 1°0. 
Converting the flux scaling factor to contrast yields 
AL’ = 10.11 + 0.25 mag for the 2021 December epoch and 
AL’ = 9.94 + 0.14 mag for the 2023 February epoch. 
Incorporating the W1 magnitude of AF Lep” 
(4.92 + 0.07 mag; Cutri et al. 2012), these correspond to 
apparent magnitudes of  L'- 15.03 + 0.26 mag and 
L' = 14.87 + 0.15 mag, respectively, and absolute magni- 
tudes of Mp = 12.89 + 0.26 mag (2021 December) and 
Mı = 12.72 + 0.15 mag (2023 February). 


3.2. TESS Light Curve 
AF Lep was observed with TESS (Ricker et al. 2015) in 
Sectors 5, 6, and 32, and processed following the procedure 


17 Here we assume that L’ — W1 = 0 since these filters are in the Rayleigh— 
Jeans tail of the F8 host star’s SED. 
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described in Bowler et al. (2023). The Science Processing 
Operations Center’s Pre-search Data Conditioning Simple 
Aperture Photometry (PDCSAP; Jenkins et al. 2016) light 
curve was downloaded using the lightkurve package 
(Lightkurve Collaboration et al. 2018). Outlier photometric 
points are removed by flattening the light curve with a high- 
pass Savitzky-Golay filter and identifying all data outside of 
three standard deviations. 

Clear modulations are evident in the final light curve with an 
amplitude of 196—296. The rotation period and uncertainty 
measurements are calculated by producing a generalized Lomb- 
Scargle periodogram (Zechmeister & Kürster 2009), fitting a 
Gaussian to the highest periodogram peak, and adopting the 
resulting mean and standard deviation. This yields a raw 
periodicity of 1.007 + 0.007 days. However, starspots on AF 
Lep could be located at mid-latitudes rather than the equator. If 
differential rotation is present, then this will impact the 
interpretation of these modulations as a rotation period. To 
account for this, we follow Bowler et al. (2023) by inflating the 
uncertainty assuming a solar-like absolute shear of 0.07 rad day !. 
Adding this uncertainty from potential differential rotation 
in quadrature with the uncertainty from the periodogram 
analysis gives a final adopted value for the rotation period 
of 1.007+0.009 days. This is similar to values of 
0.966 + 0.002 days from Messina et al. (2017), 1.03 + 0.03 days 
from Gaia Collaboration et al. (2022), and 1.01 + 0.05 day from 
Züfiga-Fernández et al. (2021).'® 


4. Results 
4.1. Common Proper Motion 


Our second epoch of relative astrometry from 2023 February 
shows no significant change in separation compared to the first 
epoch taken 13 months earlier (338+ 11 mas compared to 
342 + 8 mas in 2023), but a large counterclockwise change in 
P.A. (6228 + 1?3 compared to 72°0+ 1°0 in 2023). To test 
whether this is consistent with a stationary background source, 
perfect co-movement with the host star, or co-movement plus 
orbital motion, we compare our measurements with the 
predicted trajectories in Figure 2 using the host-star sky 
coordinates, parallax, and proper motion. At epoch 2023.090, 
the predicted separation for a stationary source is 
367 + 10 mas, and the predicted P.A. is 5674 + 170. Although 
the second epoch separation is similar, our measured P.A. of 
AF Lep b differs by about 9°2 between the two observations. 
This corresponds to a 2.00 discrepancy in separation and an 
11.00 discrepancy in P.A. We thus find that the common 
proper motion test strongly favors the source being a 
gravitationally bound companion rather than a background 
star. The apparent motion between the two epochs is therefore a 
result of orbital motion, not seen in separation (0.30) but 
detected at the 5.66 level in P.A. Assuming linear motion, the 
P.A. of AF Lep b is increasing by about 8° yr. 


4.2. Color-Magnitude Diagram and Inferred Mass 


We generate a M,H — L') color-magnitude diagram 
(CMD) of nearby substellar sources to compare against the 
L' photometry of AF Lep b using The UltracoolSheet 


18 Note that the Züfiiga-Fernández et al. (2021) measurement only used Sector 
5 of the TESS light curve of AF Lep, which may contribute to its higher 
uncertainty. 
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Figure 2. Top: test for common proper motion based on the two epochs of 
relative astrometry between AF Lep A and b. The left panels show the 
predicted relative motion of a stationary source in separation (top) and P.A. 
(bottom) with respect to the first epoch of relative astrometry (blue circles). 
Gray shaded regions represent lo and 2c uncertainties. The predicted position 
at 2023.090 is plotted as an open triangle, and the orange triangles are our 
measurements. Our 2023 observations significantly disagree with a stationary 
background source model at the 11o-level in P.A. The right panel shows the 
same comparison but in AR.A. and Adecl. instead of separation and P.A. 
Bottom: comparison of our L’ photometry of AF Lep b to the M; {H — L’) 
CMD. The gray regions show the 1c and 2c uncertainties on the L’ magnitude 
from our 2023 imaging. The background points show L (red) and T (blue) 
dwarfs from The UltracoolSheet (Best et al. 2020). Sources marked as young or 
low surface gravity are plotted in yellow. Young objects in the compilation are 
sources with independent age constraints «300 Myr. Low surface gravity 
sources have spectroscopic features of youth based on near-infrared 
spectroscopy. We also add substellar companions with H and L' photometry 
compiled within the species package (Stolker et al. 2020). The L’ 
photometry of AF Lep b implies an early T-dwarf spectral type. 


(Best et al. 2020). Single objects in the compilation with W1- 
band magnitudes, parallaxes, and H-band magnitudes in the 
Maunakea Observatories (MKO) filter system (Tokunaga & 
Vacca 2005) are selected. We convert the W1-band magnitudes 
to L'-band magnitudes using the My,—L’ relation derived in 
Franson et al. (2023). For the absolute L’ magnitude of AF Lep 
b, we adopt the value from our 2023 February imaging of 
Mz! = 12.72 + 0.15 mag. Our CMD is shown in the lower 
panel of Figure 2. The photometry of AF Lep b would imply an 
early T-dwarf spectral type although at the young age and low 
gravity of the companion, thick clouds may significantly 
redden the spectrum (e.g., Liu et al. 20162; Faherty et al. 2016). 
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From our absolute L’ photometry, we can determine the 
model-inferred mass of AF Lep b given its age of 24 + 3 Myr. 
Our approach is to draw from the age and absolute L’ 
magnitude distributions 10° times and interpolate a given hot- 
start model grid to determine the corresponding mass 
distribution. We perform this procedure for both Cond (Baraffe 
et al. 2003) and the ATMO-2020 (Phillips et al. 2020) models. 
This yields hot-start model-inferred masses of 4.53*)'3 Miup 
for Cond and 5.4 + 0.5 Mjyp for ATMO-2020. We can also 
infer masses from the Burrows et al. (1997) and Saumon & 
Marley (2008) grids by determining the bolometric luminosity 
of AF Lep b and interpolating the models. We base our 
bolometric correction on HR 8799 b, which has a similar L’ 
absolute magnitude to AF Lep b (Mj; = 12.66 + 0.11 mag; 
Marois et al. 2010). Following the procedure outlined in Brandt 
et al. (2021a), the bolometric luminosity of HR 8799 b is 
log(L/Lpo) = —4.79 + 0.10 dex, which corresponds to a 
bolometric correction BC(L/) = —0.68 + 0.27 mag. Applying 
this correction, the bolometric luminosity of AF Lep b is then 
log(L/Lyo) = —4.81 + 0.13 dex. This produces inferred 
masses of 62*17 Myy for Burrows et al. (1997) and 


5507 Myyy for the Saumon & Marley (2008) grid with the 
hybrid prescription for clouds. The cloud-free and cloudy 
versions of the Saumon & Marley (2008) models yield masses 
of 6.1483 Myup and 6313 My, respectively. 


4.3. Orbit Fit and Dynamical Mass 


Here, we perform a joint orbit fit with orvara (Brandt et al. 
2021b) of our Keck/NIRC2 relative astrometry, archival RVs, 
and the HGCA proper motions. The archival RVs were 
published in Butler et al. (2017) and consist of 20 measure- 
ments with the HIRES spectrograph (Vogt et al. 1994) on Keck 
I. The observations were taken from 2002 September to 2013 
November and have a median reported uncertainty of 
63.2 ms ! and an rms of 188.7 m s™'. 

orvara uses the parallel-tempered MCMC ensemble 
sampler in emcee to infer the orbit element posteriors. The 
following orbital elements are directly fit: companion mass 
M.omp, primary mass Myo, semimajor axis a, inclination io, 
longitude of ascending node Q, longitude at the reference 
epoch of 2010.0 Aj; and RV jitter ogy. Eccentricity e and 
argument of periastron w are parameterized as ve sinw and 
Je cosw to avoid the Lucy-Sweeney bias against circular 
orbits (Lucy & Sweeney 1971). orvara analytically margin- 
alizes over parallax, barycentric proper motion Ha and us, and 
RV instrumental zero-points. We adopt uniformative priors for 
all quantities except host-star mass, where we take a prior of 
1.2 + 0.06 Mo from its mass in Kervella et al. (2022). We use 
20 temperatures, 100 walkers, and 5 x 10? total steps (5000 
steps per walker) in the orbit fit. The first 2096 of each chain is 
discarded as burn-in. Convergence is assessed by verifying that 
multiple runs with different starting parameters produce 
consistent parameter posteriors. 

Table 2 shows the results of our orbit fit. Figure 3 displays 
the posterior distributions for selected orbit elements alongside 
a comparison of the relative astrometry of AF Lep b to a swarm 
of orbits drawn from the fit. We measure a dynamical mass of 
3974 Myup, semimajor axis of 8.4*11 au, orbital inclination of 
50*75°, and eccentricity of 0.24*077. The orbital period is 
22 4: 5 yr, and the time of periastron is T) = 2014.1*2J. AF 
Lep b and its host star have a mass ratio of 0.0026 + 0.0005. 
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Table 1 
Properties of AF Lep A and b 
Property Value References 
AF Lep A 
O/000.0 05:27:04.76 1 
62000.0 -11 :54:03.5 1 
la? (mas yr~') 16.915 + 0.018 1 
ls (mas yr~!) —49.318 + 0.016 1 
a? (mas yr ?) 0.0191 + 0.0024 3 
a? (ms ‘yr )) 2.43 + 0.30 3 
7 (mas) 37.254 + 0.020 1 
Distance (pc) 26.843 + 0.014 1 
SpT F8V 2 
Mass (Mo) 1.20 + 0.06 3 
Radius (Ro) 1.25 + 0.06 4 
Age (Myr) 24 t3 5 
Tere (K) 6130 + 60 6 
[Fe/H] (dex) +0.19+ 0.02 7 
vsini (km/s) 50+5 8 
log(R'uk) (dex) -4.27 9 
log(Lx/Lpo1) (dex) —3.48 10 
Prot (d) 1.007 + 0.009 11 
ix (°) 54)! 11 
RUWEpr3 0.918 1 
B (mag) 6.83 + 0.02 12 
V (mag) 6.30 + 0.01 12 
Gaia G (mag) 6.210 + 0.003 1 
J (mag) 5.27 + 0.03 13 
H (mag) 5.09 + 0.03 13 
Ks (mag) 4.93 + 0.02 13 
W1 (mag) 4.92 + 0.07 14 
AF Lep b 
2021.970 (mas) 338 + 11 11 
42021.970 C) 62.8 + 1.3 11 
P2023.090 (mas) 342 +8 11 
2023.090 C) 72.0 + 1.0 11 
AL’ (mag) 9.94 + 0.14 11 
L’ (mag) 14.87 + 0.15 11 
Mz, (mag) 12.72 + 0.15 11 
log(L/Lyg) (dex) —4.81 + 0.13 11 
Dynamical Mass (Myup) 29907 11 
Semimajor Axis (au) 8.4414 11 
e 0.247071 11 
i, (°) 50*2, 11 
Period (yr) 22+5 11 


References. (1) Gaia Collaboration (2022), (2) Gray et al. (2006), (3) Kervella 
et al. (2022), (4) Kervella et al. (2004), (5) Bell et al. (2015), (6) Ammons et al. 
(2006), (7) Gaspar et al. (2016), (8) Glebocki & Gnacinski (2005), (9) Isaacson 
& Fischer (2010), (10) Wright et al. (2011), (11) this work, (12) Hgg et al. 
(2000), (13) Skrutskie et al. (2006), (14) Cutri et al. (2012). 

è Proper motion in R.A. includes a factor of cos 6. 

> Calculated from the proper motion difference between the Hipparcos-Gaia 
joint proper motion and the Gaia EDR3 proper motion in Brandt (2021). 


To assess the potential of systematics in the relative 
astrometry or Hipparcos-Gaia proper motions impacting the 
companion's dynamical mass, we perform a suite of additional 
orbit fits to synthetic data with systematic offsets from the 
nominal measurements at the 2c and 4c levels. Each of these 
fits uses 20 temperatures, 100 walkers, and 10° steps. We first 
vary each epoch of the relative astrometry by +20 and +40 
while keeping the other epoch constant. Here, the separation 
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Table 2 

AF Lep b Orbit Fit Results 

Parameter Median +10 95.4% C.I. Prior 
Fitted Parameters 
Meomp (Map) 3.2107 (2.1, 4.7) 1/Meomp (log-flat), Mcomp € (0, 10? MoJ? 
Mros (Mo) 1.20 0.06 (1.08, 1.32) 1.20 + 0.06 M (Gaussian) 
a (au) 8.4414 (6.2, 12.8) 1/a (log-flat), a € (0, 2 x 10° au]? 
i, ©) 50+% (23, 66) sin(i), 0° < i < 180° 
Je sinw —0.1*02 (—0.8, 0.7) Uniform, — Ve < Je sinw € Je 
Je cos w —0.16*050 (—0.60, 0.58) Uniform, — Ve € Je cosw < Je 
Q (9) 66+13° (10, 122)" Uniform, —180? < w < 540° 
dre CY 160720 (10, 350) Uniform, —180? < Ayer < 540? 
Parallax (mas) 37.254 + 0.020 (37.214, 37.293) 37.254 + 0.020 mas (Gaussian) 
Ha (mas yr) 17.111 + 0.024 (17.074, 17.154) Uniform, Ha € (—oo , oo )* 
us (mas yr~') —49.183*0021 (—49.211, —49.142) Uniform, jus € (—00 , oo Y 
RV Jitter ogy (ms) 17733 (124, 265) 1/opy (log-flat), ory € (0, 10° ms ^! 
Derived Parameters 

P (yr) 225 (14, 42) 
e 0.247027 (0.02, 0.71) 
w (°) 210798, (20, 340) 
Ty (JD) 24567007898, (2448400, 2458500) 
To (yr) 2014.1*24 (1991.3, 2019.1) 


0.0026 + 0.0005 


q (= comp/Mhost) 


Notes. 


(0.0017, 0.0037) 


* The posterior distribution for © consists of two distinct peaks separated by 180°. The values shown in the table correspond to the lower peak. The other peak is 


located at 250*$0^ with a 95.4% confidence interval of (170, 350). 
> Mean longitude at the reference epoch of 2010.0. 
* Note that these are improper priors since they do not normalize to one. 


and position angle are both offset by the same amount for each 
instance (for example, +40 in separation and +4c in position 
angle). We find that the dynamical mass is resilient to 
systematics at the 2c and 4c levels for both epochs, with all 
resultant companion masses being consistent with our dyna- 
mical mass to within 1c. The companion mass posteriors from 
the mock runs have similar uncertainties to the dynamical 
mass. Changes to the HGCA proper motions have a larger 
impact on the dynamical mass. Shifting the Hipparcos-Gaia 
proper motion by —2c or the Gaia proper motion by +20 
decreases the dynamical mass by about 1c from our nominal 
case to ~2.5 Myyy. Applying the opposite offset increases the 
dynamical mass by about lo to ~3.9 Mjup. Offsets at the 4o 
level change the dynamical mass by about 2c to about 1.7 Myup 
or 4.5Mjyp. We thus find that the dynamical mass is more 
significantly affected by systematics in the Hipparcos-Gaia 
proper motions than the relative astrometry. 

Our dynamical mass measurement is slightly lower than 
model-inferred masses from the planet’s L’ photometry and 
associated bolometric luminosity. It is 1.70 lower than the 
Cond-inferred mass of 4.53*032 Myup, 2.60 lower than the mass 
inferred via ATMO-2020 of 5.4 + 0.5 Myyp, 2.70 lower than 
the mass inferred from Burrows et al. (1997) of 6.2799 My, 
and 2.40 lower than the Saumon & Marley (2008)-inferred 
mass with the hybrid prescription for clouds. Here, one-sided 
Gaussian-equivalent c values are determined following Fran- 
son et al. (2023). All inferred masses are generated using hot- 
start models, which treat planets as having arbitrarily large radii 
and high entropies at early ages as an initial condition. Cold- 
(e.g., Marley et al. 2007; Fortney et al. 2008) and warm-start 
(e.g., Spiegel & Burrows 2012; Marleau & Cumming 2014) 


models start with lower amounts of initial entropy to emulate 
the loss of energy through accretion shocks in planet formation. 
Since this lowers the luminosity of planets at young ages 
relative to hot-start models, the discrepancy between our 
dynamical mass and model predictions is larger if a cold- or 
warm-start model is assumed. 

One way this modest discrepancy can be reconciled is if the 
planet is younger than the star, perhaps due to delayed 
formation in a disk. This would cause the planet to be brighter 
than expected based on the host-star age. For the case of AF 
Lep b, the planet would require an age of 16*7 Myr for the 
Cond mass to match the dynamical mass. Other models require 
lower masses; for the prediction from the hybrid Saumon & 
Marley (2008) model to match the dynamical mass, the planet 
would need to be 9*3 Myr. This would imply a formation 
timescale of ~5-15 Myr, which would suggest formation via 
core accretion, which takes many Myr to operate, over rapid 
giant planet formation via disk instability. Typical protoplane- 
tary disk lifetimes are 5-10 Myr for Sun-like and lower-mass 
(<2 Mọ) stars (Pfalzner et al. 2022). Other possibilities for 
resolving this discrepancy include a younger age of the system, 
young L/T transition objects being more luminous than 
predicted by current evolutionary models, and another planet 
in the system biasing the dynamical mass measurement. 


4.4. Planet—Disk Interactions 


There are several potential avenues for an inner planet to 
dynamically interact with an outer debris disk. One possibility 
is the planet truncating the disk's inner edge. Due to the close 
separation of AF Lep b from its host star, it is unlikely that the 
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Figure 3. Joint posterior distributions of Mnost, Mcomp» 4, e, and i for the orbit fit of AF Lep b. The diagonal panels show the marginalized distribution for each 
parameter. The contours on off-diagonal panels denote 1c, 20, 3c, and 4c levels in the joint posteriors. The upper right plots compare the separation and position 
angle of AF Lep b over time against a swarm of orbits drawn from the orbit fit, shown in gray. The black curves highlight the maximum-likelihood orbit. 


planet is truncating the disk from its current orbit. Pearce et al. 
(2022) found that the disk inner edge could be truncated by a 
1.1 + 0.2 Myyp planet at a semimajor axis of 35 + 6 au. Larger 
objects at closer-in semimajor axes could also produce the 
same effect, but the dependence between mass and semimajor 
axis is sufficiently steep that by the 8 au semimajor axis of AF 
Lep b, only stellar companions are capable of truncating the 
disk (see Equation (6) of Pearce et al. 2022). This remains the 
case if we allow for the possible eccentricity of AF Lep b and 
the uncertain location of the unresolved disk. If the debris disk 
was truncated by planets, then either additional, unseen planets 
are responsible or AF Lep b was closer to the disk in the past 
and has since migrated inwards. 


While AF Lep b is unlikely to be truncating the disk, its 
possible eccentricity could drive the disk into an asymmetric 
shape. Equation (5) in Pearce et al. 2022 relates the semimajor 
axis and eccentricity of a planet's orbit to the eccentricities of a 
debris disk's edges. A planet at 8au with an eccentricity of 
0.25 would drive a narrow disk at 46 au to an eccentricity of 
~0.1. The disk is currently unresolved, so no eccentricity or 
width information is available. Future resolved observations of 
the disk have the potential to determine the disk's eccentricity 
and directly examine the interaction between an eccentric 
planet and its debris disk. 

A third possibility is stirring, where the planet excites 
planetesimals in the disk, causing them to collide and release 
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the observed dust. Equation (15) of Mustill & Wyatt (2009) 
gives the mass of an interior planet required to stir debris as a 
function of planet eccentricity and location. For an eccentricity 
of 0.25 and semimajor axis of 8 au, planets with masses greater 
than 0.45 Mj, are able to stir the AF Lep debris disk. AF Lep 
b, with its mass of 3.210? Myyp, may be actively stirring the 
disk. This is significant as there is growing evidence that many 
debris disks are unable to stir themselves via self-gravity 
(Krivov & Booth 2018; Pearce et al. 2022, though see also 
Najita et al. 2022). Planetary stirring could therefore be 
prevalent in generating the observed dust in debris disks. Note 
that if AF Lep b has a small but nonzero eccentricity, it could 
still stir the disk, provided the eccentricity is above ~0.04. This 
threshold is determined by rearranging Equation (15) of Mustill 
& Wyatt (2009) to give the minimum eccentricity capable of 
stirring the disk as a function of planet mass. 


4.5. Stellar Inclination and Minimum Obliquity 


The line-of-sight inclination of AF Lep, ią, can be inferred 
from its rotation period, stellar radius, and projected rotational 
velocity. We follow the Bayesian framework of Masuda & 
Winn (2020) to derive ią using the analytical expression for its 
posterior distribution from Bowler et al. (2023, their Equation 
(9)), which accounts for the correlation between projected and 
equatorial rotational velocities. Here we adopt the TESS 
rotation period of 1.007 + 0.009 d, the v sini, value of 50 + 
5 km s! from Glebocki & Gnacinski (2005), and radius of 
1.25 +0.06 Ro from Kervella et al. (2004). The maximum 
a posteriori value is 54°. The 68% credible interval spans 45?— 
65? and the 95% credible interval spans 39?-82^. 

The stellar inclination together with the inclination of the 
orbital plane of AF Lep b (i,) provides information about the 
minimum stellar obliquity—whether the spin axis of AF Lep is 
aligned or misaligned with respect to the planet. Our orbit fit 
indicates an inclination of i, = 50*?,°. This is in excellent 
agreement with the stellar inclination of ią = 54*5'°. The polar 
position angle of AF Lep—which describes its absolute 
orientation in the sky plane—is not known and can therefore 
not be compared with the longitude of ascending node of the 
orbit to assess the true obliquity angle. However, the agreement 
of the inclination angles is nevertheless consistent with spin- 
orbit alignment, which reinforces trends seen in other systems 
with imaged planets but differs from the more diverse 
architectures of brown dwarf companions (Bowler et al. 
2023). The debris disk around AF Lep has not been resolved; 
if future observations resolve the disk, its inclination could be 
compared with the stellar and orbital inclinations to further 
assess angular momentum alignment in the AF Lep system. 


5. Summary 


AF Lep b is a young giant planet directly imaged inside the 
debris disk of a Sun-like star in the 24+ 3 Myr ( Pic moving 
group. The system was identified as a promising target based 
on the low-amplitude (2.43 20.30 ms ! yr )) astrometric 
acceleration the planet induces on its host star. Our main 
conclusions are summarized below: 


1. Our two epochs of NIRC2 observations in L’ with the 
VVC separated by 13 months confidently show that AF 
Lep b is gravitationally bound and is undergoing orbital 
motion around its host star. We can rule out AF Lep b 
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being a background star at the Ilo level, and orbital 
motion is detected at the 6c level. 

2. AF Lep b has a contrast of AL’ = 9.94 + 0.14 mag at 
0734 and an absolute magnitude of Mr, = 12.72 + 0.15 mag. 
Based on the age of the host star, and assuming the planet 
was formed at the same time, the implied planet mass from 
hot-start evolutionary models is 224-6 Mjyp. 

3. Our orbit fit incorporating the relative astrometry, 
astrometric acceleration, and RVs yields a dynamical 
mass of 3.20; Myup, a semimajor axis of 8.4*]:3 au, and a 
modest eccentricity of 0.24*077 for AF Lep b. The mass 
of AF Lep b is too low at an orbital distance of 8 au to 
truncate the outer debris disk, which is located at 
7-46 + 9 au (Pearce et al. 2022). Another massive planet 
or series of low-mass planets may be present at wider 
separations. 

4. The dynamical mass (3.204 Myup) is somewhat lower 
than predicted masses (%4—6Myyp) from hot-start 
models assuming a nominal age of 24 + 3 Myr from the 
host star. This could be caused by a slightly younger age 
for the host star, delayed formation in a disk, systematic 
uncertainties in the evolutionary models, or another 
companion in the system. 

5. The stellar inclination (ią = 5414!°) and orbital inclination 
(i, = 50725?) are in good agreement, which is consistent 
with the system having spin-orbit alignment although 
polar position angle of the star is unknown, so this 
agreement represents a lower limit on the true obliquity. 


AF Lep joins other young planet hosts with debris disks such 
as B Pic, HR 8799, HD 206893, and HD 95086, reinforcing 
indications of a higher frequency of long-period planets 
orbiting stars hosting debris disks (Meshkat et al. 2017). 
Within the 8 Pic moving group, the mass of AF Lep b is higher 
than 51 Eri b (z:2 Mjup; Macintosh et al. 2015) and lower than 
p Pic b and c (8-9 Mj,y; Nowak et al. 2020), contributing to a 
growing sequence of giant planets with common ages. AF Lep 
b is the lowest-mass imaged planet with a dynamical mass. We 
expect it to become an excellent target for follow-up orbital and 
atmospheric characterization. 
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